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Square Root Auxiliary for Routine A7, Floating Decimal 

Arithmetic Routine (DOI or SADOl) 

Closed 

30 

S3 ^ Location of Floating Accumulator 

S^, Location of first word of Routine A7 

0, 1, 2. all of Floating Accumulator 

Approximately 25 msec 

This routine computes the square root of the number in 

the floating accumulator to full double precision, using 

Newton's method,, as described in Routine Rl - 116. The 

contents of the Floating Accumulator are standardized 

by this routine on entry , so the programmer need not use 

an interpretive N2F instruction for this purpose before 

entry » The routine is entered with a 8j x F order, where 

x is the location of the first word of this routine , and 

control is returned to A'7 so that the next interpretive 

instruction obeyed will be the one following the 8j order, 

with the computed square root in the floating accumulator „ 

If the floating accumulator is negative , this routine 

will stop on a division hangup. 

If AT is operating in the fixed-point mode, an additional 

right shift of the floating accumulator of 256 - N;(2S3) 

places is required for an arithmetically correct result. 
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00K(RA2) 









L5 3L 








^2 25L 


Set shift address 




1 


L5 S3 
^0 3S3 






2 


50 1S3 








F5 2L 


Standardize Floating Accumulator 
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26 284S4 
10 IF 






k 


kO S3 
L5 2S3 
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LO 5S3 








Lh 387S4 


Correct exponent for standardization 
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10 IF 








kO 2S3 


and set bias for new exponent 
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S5 F 








00 IF 


Test for odd exponent 
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36 11L 








F5 2S3 


jump if was even 
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kO 2S3 








L5 S3 


Correct exponent and Floating Accumulator 
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10 IF 
kO S3 






11 


S5 F 
kO 1S3 






12 


51 S3 








10 IF 


Compute single-length square root 
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SJ F 








kO 2F 


y Q = (1/2 + l/2x) 
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L5 S3 
50 1S3 






15 


66 2F 






1 


S5 2^L 







LOCATION 
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ORDER 



LO 2F 
32 18L 
10 IF 
Ik 2F 
22 13L 
FS 8F 
36 23L 
F5 25L 
42 25L 
L5 2F 
40 3S3 
50 386S4 
L5 15L 
26 244S4 
L5 2F 
80 IF 
26 21L 
L5 S3 
50 1S3 
10 IF 
00 IF 
40 S3 
S5 F 
40 1S3 
L5 2F 
50 386S4 
10 IF 
26 n6s4 



NOTTS. 



PA.GE .2 



RA. 2 



n /rt /Floating Accumulator N 

y^n = y v + !/ 2 (- — — ^r — +; -y k ) 



k+1 "k 



y* 



Test ly, I > 1/2? jump if not, 



set shift to scale down 



double precision division to form 



F. A, 



n 



Scale up divisor 



properly scale 



Floating Accumulator 



n 



TJl A 

re-enter A7 to form l/2 (^ * > *■ '+ y ) 
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